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ABSTRACT 



Three independent wave celerity data sets, measured on natu- 
ral teaches, are compared with linear, bore, solitary, and 
hyperbolic wave theories. In the range or relative water 
depths (.0C6 < h/T Z < 13 cn/s~) and wave heights (.1 < H/T~< 
3 cm/s ) tested, hyperbolic wave theory, which is an assymp- 
totic form cf cncidal theory in shallow water, agreed most 
closely with measured wavs celerities. Linear wave theory 
also gave satisfactory results; but bore and solitary wave 
theories overestimated the observed wave speeds. It is con- 
cluded that the observed waves are weakly dispersive in 
amplitude and that care must b° taken to apply the theories 
only in their regime of validity. 
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I. INTRODUCTION 

Considering all the various theories available to 
describe wave celerity, surprisingly little work has been 
oriented toward comparing wave theories with data, particu- 
larly field data. Wave celerities acquired in the labora- 
tory were tested against wave theories by LeMehaut* et a 1 
(1968); they concluded that cnoidal wave theory as proposed 
by Keulegan and Patterson (1940) , with equations and tables 
by Masch and Wiegel (1961), gave a 'best fit* solution. 
Other laboratory studies have suggested appropriate wave 
celerities to be described by solitary wave theory (Ippen 
and Kulin, 1955) and (Kishi and Saeki, 1967), Stokes theory 
(De, 1955), by per tolic/cn cida 1 theory (Iwagaki, 1968) . Cel- 
erities measured in the field have been compared with linear 
theory (Thornton and Guza, 1982) , and bore theory (Bradshaw, 
1982) and (Suhayda and Pettigrew, 1977) . The validity of 
various theories depends upon the relative depth (h/L) , 
where h is the water depth and L is the wavelength, arc. the 
relative waveheicht (H) measured in terms of a wavelength 
(H/L) or depth (H/h) . 

In this study, various celerity formulae given by pro- 
gressive wave theories are tested against celerities 
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measured in the field. To test the theories against data, 
the phase speed equations require seme combination of water 
depth, wavsheight, and wave period (T) . Several relatively 
large field data sets did not provide one or the other cf 
these parameters and could not be used. Data which repre- 
sents reflected waves from a sceep- beach such as a data set 
collected at Fort Ord Beach, California (Sallenger et al, 
1983) could net be used. Also, another large data set, col- 
lected by University of California, Berkeley, California 
(Hoffitt, 1953) using photographic methods, was no": included 
in the study. This is because the waveheight was omitted in 
the measurements. Waveheight must be included tc test 
higher order theories, especially if amplitude dispersion is 
to be examined. 

Five wave theories were compared with wave celerities 
measured in the field in water depths from .07 to 10.0 
meters, fer wave heights from .12 to 1.0 meter, and for 
periods varying between 4,7 and 18.3 seconds. Stokes (third 
order) theory as presented by Hunt (1953) proved unsuitable 
due tc the relatively shallow water depths encountered. The 
four remaining theories tested are: linear theory; solitary 
theory as modified by Laitone (1959) ; bore theory as given 
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by Keller et al (196 0), and hyperbolic theory of iwagski 
(1968), which is an assymptotic form for cnoiaal theory in 
shallow water. 

The data have been plotted as a function of waveheight 

and water depth versus period in Figure 1. This chart, 

adapted from LeMehaute (-1 976), attempts to quantify regimes 

of applicability. For example, the diagram indicates that 

the data dees not fall in the Stokes (third order) theory 

regime where experimentation proved to be true. Additicn- 

2 3 

ally, the Ursell (Ur) parameter (Ur = HL /h ) , plotted as a 
dashed line, is commonly used to parameterize the nonlinear 
waves in shallow water (Thornton and Guza, 1982). For 
Stokes theory to be valid, the Ur se 11 parameter should be 
small, which is not the case for most of the data consid- 
ered. In the case of very long waves in shallow water, the 
Ursell parameter becomes meaningless since it is directly 
proportional to the square of the wavelength. 

To determine the range of validity for various theories, 
the underlying assumptions are examined. All the theories 
presented assume: the motion to be irro rational , the fluid 

to be incompressible , no mean current flow, and normal wave 
incidence. As a consequence of the wave theory assumptions. 
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the Ursell parameter can now be quantified and is ir.c 
for comparison purposes. The assumptions for the va 
theories and their celerity equations are as fellows: 



LINE J F : 



c - (gk tanh kh) 1/2 

Assumes: H/h << 1, H/I << 1, h/L << 1, and Ur >> 1 for 
low water, h/I >> 1 and Ur « 1 for deep water. 



SOLITARY 



/gh fl.O + I (§) - 2 — (^-) 2 1 
[ 2 v h' 20 v h' J 



Assumes: H/h < 1, h/L < .1, Ur=0(1). 



90 RE: 

c - /jr ci.o + f) l/2 d.o + f^, 1 ' 2 

Assumes: H/h=C(1), h/L << 1, Ur >> 1 (very shallow wate 



1 u d e d 
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HY fE R EOLI C : 



c 



/ gh ~ ( l.o 



( 1.0 



- 2th 5 [»-° + <x -° + fh> 
+ 2H) (II ) 2 I I (I - I) - 

+ Kh ; | K ' K 4 ' 





Jssuiss: E/h < 1, h/L < .1, Ur >> 1, where K is d~ 

below. Iwagaki (1968) linearized rhe computationally c 
cult Jacobian, elliptic function ir.ro a modulus K , when-; 

K . T/gTC ST ( H)l/2 r i0 . a ( H >n Jm 

in which a=1.3, n=2 and m=.5 for H/h - 0.55 
and a=0.5u, n=1.5 and m=1 for H/h > 0.55 
This approximation greatly simplified the analysis. 



<“> 



fined 



(5) 
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II 



. EXPERIMENTS 

Three field data se:s are used to compare with various 
theories. Experiments conducted at Torrey Pines B^ach and 
Leadfcetter Beach, California were both part of the Nearshore 
Sediment Transport Study (NSTS) (Seymour and Duane, 1 976). 
These beaches were chosen for their relatively simple beach 
plan, both unbarred, and essentially straight and parallel 
nearshore contours. The third data set was from Seven Mile 
Beach, Australia. 

A., SEVEN KILE EEACH, AOSIBALIA 

The Sever. Mile Beach (Shoal haven Bight) experiment on 
the south coast of New South Wales, Australia, was conducted 
in early 1982 by Eraashaw (1982) . He examined the relation- 
ship between bore velocity, height, and water depth by ana- 
lyzing movie camera pictures using stakes driven into the 
sand as references. This technique examines each individual 
wave by counting the number of m c vie frames between the four 
reference stakes to determine wave speed while reading bcr= 
height and water depth directly from the incrementally 
marked stakes. 
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The beach is composed of fine quartz sand and has multi- 
ple cffshcre tars. The beach slcpe outside was 0.04, while, 
the inner surf zone had a slope of 0.03. The waves gener- 
ally broke cn the outer tars and then reformed and propa- 
gated as teres inside the surf zone. The outside breakers 
ranged from 1.0 to 1.5 meters. 

The data are limited tc cnly shallow water bores inside 
the breaker line, where tore heights were from 0.14 tc 0.30 
meters in water depths cf 0.07 to 0.42 maters. Data points 
are plotted as circles in Figure 1 and generally reside in 
the very shallcw water, small waveheight regime. Cf the 27 
data points available, three points were discarded because 
bore height greatly exceeded the water depth (H/h > 3) . 

Zxaci per icd/f reguency data was not available for the 27 
runs; however, during the two-day experiment, an 3 to 12 
second period was observed. Hence the mean (10 second) 

peried is chcssn for computational purposes. This assump- 
tion is adequate since only shallcw water bores were consid- 
ered and they are essentially no ndispersiv e. 

B. TCRR3Y FINES EZACE, CALIFORNIA 

Field measurements were made at Torrey Pines Beach near 
San Diego, California in August 1973. A detailed 
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description cf the experiment can be found in Guza and 
Thornton (1980), Guza and Thornton (1982), and Thornton and 
Guza (1982). Tcrrey Pines is a gentle sloping (.02) beach 
and is composed cf fine grain sand. The waves were gener- 
ally narrow banded and approached the shora at a near normal 
angle. Directional properties of the waves were measured 

using a linear, five pressure sensor array in 10 meters 
depth. The angle of swell approach was limited to the maxi- 
mum and minimum cf +/- 15 degrees due to sheltering by off- 
shore islands and coastline restrictions, but the angle was 
generally less than +/- 5 degrees. Refraction analysis 
shewed that the predominant swell waves (T = 13 sec.) start- 
ing at 15 decrees in 10 meters depth results in angles cf 
incidence cf 6.5 degrees in 3 meters depth and 4.9 degrees 
in 1 meter depth (Thornton and. Guza, 1982). 

The 92 data points, plotted as +'s in Figure 1, are gen- 
erally restricted to the shallow water regime. Waveheights 
ranged from small to 2.0 meters, water depths at sensor 
locations from 0.27 tc 7.C meters, and mean periods from 9.0 
to 18.0 seconds. The experimental domain included spilling 
cr mixed spilling and plunging breakers. 
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C. LEADBETTEH BEACH, CALIFORNIA 

The Leaabetter Beach, Santa Barbara, California, experi- 
ment was conducted during the period of 30 January to 23 
February 198 0. Experimental details are described ir. Gable 

(1961). leadbstter is a relatively straight, steeper slop- 
ing (.05) beach composed of fine to medium, well-sorted 

sand. A series of storms, resulted in abnormally large 

waves and bsach erosion for the period commencing 4 Febru- 
ary, 1930. The storms formed what has been described as a 
•50 year' storm event. Measured wave heights during the 
experiment ranged from 0.18 to 1.9 meters, with sensor 
depths up to 10 meters, and periods from 4.7 to 16 seconds. 
Breaking wave types were of both the spilling and plunging 
variety. The 83 Santa Barbara data points, plotted as tri- 
angles in Figure 1, exhibit the largest waveheights and 



water depths examined. 

Leadbetter Beach has an east-west orientation which is 
counter to the north-south orientation of the California 
coastline. Ihe predominant northwest ocean swell entering 
the narrow gap between Joint Conception and the Channel 
Islands would have to refract nearly 90 degrees to approach 
normal to the beach. As a consequence, the ocean swell 
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waves approach at a relatively large, well focused angle 
from the west. At other times, storm waves generated inside 
the Channel Islands approached at large angles from the 
east. Therefore, the wave angularity has to be considered 
in the wave celerity calculations. 

Incident offshore sea-swell was measured using a fcur 
pressure sensor square array with 6 meter legs located ir. a 
water depth of 8 meters. From these sensors, the mean inci- 
dence angle at the teak frequency was determined for the 
nine data sets. The wave rays were then manually refracted, 
using Snell's Law, from the pressure array location to a 
point where the bottom contours could be considered straight 
and and parallel to the beach. Shoreward of this point, the 
incident wave angles were calculated using a constant 
refractive coefficient. Deep water angles varied between 
+/- 20 degrees and were refracted into 4 meters of water 
frcm there into the shallower water region of the current 
sensors. For example, wave angles of 20 degrees nr. 8 meters 
depth, period cf 12.0 seccnds, resulted in refracted angles 
cf 10 degrees in 3 meters decreasing to 5 degrees in 1 teeter 
of water. 
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Spectral methods vers utilized to determine mean celeri- 
ties and wave heights using current and pressure sensor data 
from San Eiegc and Santa Barbara. The celerity and wave 
height calculations are a r. average over many waves. Since 
the record lengths were 34 minutes and the mean period (T) 
was about 12 seconds, approximately 170 waves are averaged 
for the spectral estimate cf wave celerity. This is in con- 
trast to the camera methods of Bradshaw (1982) where each 
individual wave was photographed as it passed the reference 
stakes and waveheight and water depth were read directly. 

Pressure and current meter data were telemetered to 
shore where they were digitally recorded at a rate of 64 
sa mples/seccnd . The data were averaged to 2 samplss/seccnd 
which results in a Nycuist frequency of 1.0 Hertz. Records 
were then compiled into 4C96 da ta points. By breaking the 
series into 256 point records, the phase, kinetic energy, 
and coherency spectra were computed with 32 degrees cf 
freedom. 

Celerity spectra were calculated from the phase spectra 
measured between adjacent pairs cf current meters located in 
a line normal to the beach. The actual celerity (C x ) was 
computed using: 



c x (f> 



2tt f Ax 

<P (f ) 



( 6 ) 
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where, (ax) is the distance between sensors, ($) is the 
phase difference between sensors, and (f) is the frequency. 

It was shewn by Thornton and Guza (1982) , for the Torrey 
Pines data, that the waves can be considered frequency nen- 
dispersive such that the celerity spectra is constant, at 
least across the energetic region of the spectrum. Hence, a 
mean celerity, representative of the entire spectrum, was 
chosen at that value corresponding to the peak frequency in 
the energy spectrum and is used in the data comparison here. 
Also, the frequency at the spectral peak generally coincided 
with the maximum coherency. 

Hence, the Tcrrey Pines Beach celerity data were calcu- 
lated using (6), where f = f(peak). At Santa Barbara the 
wave angularity was taken into account. Since the celerity 
spectrum was calculated using instruments in a line normal 
to the beach, only the x component of celerity is measured 
directly. The total ’mean' phase speed is calculated using: 





(7) 



cos a (f ) 
P 




frequency <f p ) - calculated as explained amove. 
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The root mean square cf 



waveheights was calculated from 



the sea surface variance t! 0 , where: 



H RMS ~ 1/8 M o 



(3) 



The sea surface 


elevation 


variance is 


proportional to the 


potential energy 


« 


It is 


assumed that 


the potential and 


kinetic energies 


are 


equal 


in the wave 


system. Hence, K 0 


was calculated 


using 


the 


current meter 


kinetic energy cr 



using the potential energy of the pressure sensors by apply- 
ing the appropriate linear theory spectral transfer func- 
tion . It is emphasised that when utilizing spectral 
decomposition cf the current sensors, the data was averaged 
in time such that the waveheights are representative ov-r 
the entire time cf the measurement. 

Mean water depths were measured from pressure sensors 
associated with the current meters. 

The error associated with Tcrrey Pines Beach celerity 
calculations is estimated to be within 5% for the deep water 
sensors and decreases to 1.2% shoreward (Thornton and Guza, 
1982) . The error estimates are based on the stability cf 
the phase estimates and measurement errors. Estimating sen- 
sor spacing errors, depth errors, and refraction problems 
account for the larger error value in deeper water. 
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Additional errors can fee attributed to the modification 
of the phase speed due to mean currents. The mean currents 
were net subtracted out during spectra computations, hence x 
(off shore) ar.d y (longshore) currents must be addressed. 
Lsadtetter Beach, with the largest angles of swell inci- 
dence, exhibited a nominal longshore component of 40 cm/sec. 
(max cf 80 cm/sec) in 1 meter of water, decreasing seaward, 
to 8 cm/sec. in 2 meters cf water. Since maximum incidence 
angles were or. the order of 10 degrees in 3 meters ar.d 5 
degrees in 1 meter, the error associated with the measured 
celerity is on the order cf 0.5% (1.0% maximum) and hence it 

is negligible. The current meters were located in the lower 
half cf the water column. The measured mean currents gener- 
ally exhibited a slight net offshore flow which is presumed 
to be balanced by the onshore mass transport in the crest- 
trough region. Therefore, the net on-offshore mass trans- 
port is expected to be significantly small and its effect on 
i he celerity can be considered negligible (Thornton and 
Gu za , 1982). 
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RESULTS 



The three data sets are compared in Figure 1 showing the 
approximate regions of validity for the various theories. 
The data do net compare well with the breaking criterion for 
solitary waves of H/h = C.78. Note that approximately 15 
data points fall outside the H/h=0.78 curve. These particu- 
lar waves were part of a study looking at bores during their 
run-up excursion (Bradshaw, 1982). The waves have exceeded 
the maximum limit for a solitary wave, but bores do net nec- 
essarily compare with solitary wave theory. Nevertheless, 
studies have indicated that bores should be classified 
according tc their height-tc- depth ratio (H/h) . Peregrine 
(1966) places beres into three categories. For H/h<0.28, 
the bore is called undular and is composed of a series cf 
undulations radiating behind the leading wave. For 
0. 28<H/h< 0 .75 , bores are considered partially developed and 
may break, while, the bore is termed fully developed and 
subject tc large-scale turbulent breaxing for H/h>0.75 
(Suhayda and Pettigrew, 1977) . Conceivably, Figure 1 should 
be modified in the very shallow water regime tc account for 
the three bore categories. 
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The data for Torrey Fines and Leadbezter Beaches also 
depart from the H/h = 0.76 slope. At Torrey Pines Beach, 
Thornton and Guza (1982) found waves inside the surf zone an 
saturation, H (EM S) / h = 0.44. A possible difference is tha- 
tha HiP.MS) waveheight, calculated for Torrey Pines and Lead- 
bezter Beach data do not necessarily correspond to the indi- 
vidual waveheight used in :he solitary wave criteria. 
Individual waveheight is also used for Seven Mile Beach 
data . 

The measured phase speeds from the three experiments are 



compared 


with the four wave 


Jn 0 


cries in 


figures 2-5 


solid lin 


e in figures 2 thru 


5 


denotes 


a perfect f 



(slope=1.0), the dashed is the ’best fit' linear equation 
line, and the dotted lines are the 95% confidence intervals. 
Table lisa synopsis of a linear, least squares, regression 
comparing measured values versus the theoretical celerity. 
The equations cf 'best fit' lines are given as well as the 
correlation coefficient (Corr). All data points are unbi- 
ased (given equal weight) when computing Table 1 values. 

The 953 confidence intervals (dotted lines) have been 
placed on Figures 2-5 for statistical purposes. Any pre- 
diction about an individual Y (actual) associated with a 
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given X (theory) will be most meaningful near the mean of X. 
Hence , the 95?? confidence limits are bowed inward (toward 
the regression line) near the -mean of X. The curvature 
associated with the confidence lines is slight due to the 
significant number (199) cf data points. 

A perfect fit of the data to a theory consists of a 
regression line whose slope is 1.0 and intercept is 0.0. 
The regression equation which is the best approximation cf 
the perfect fit is presumed to be the most accurate. It is 
apparent from the 'best fit' (dashed) lines of figures 2-5 
that hyperbolic wave theory most closely approximates the 
perfect solution (solid line) , while linear wave theory is 
only slightly less accurate. Bore theory appears to give a 
near constant c ve r-e s timaticn of phase speed whil e solitary 
theory does well at the lew end (<500 cm /sec) cf the spec- 
trum. Statistically, the errors tend to nullify each ether 
when the data sets are combined. This is a recognized prop- 
erty cf a simple linear regression and becomes quite evident 
as the number cf data points (N) approaches infinity. 

While the equations of the least squares linear regres- 
sion suggest that hyperbolic wave theory may have a slight 
edge; the same cannot be said if each beach is analyzed 
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separately. Torrey tines Beach is best explained by here 
theory. Seven Mile Beach ty solitary theory, ar.d Leadbetter 
Beach by hyperbolic theory. We also notice at Seven Mile 
Be ac h that nene cf the theories did particularly well at 
predicting wave speeds. This should be no surprise since 15 
cf the 2 u data points fell outside the regime cf the wave 
theories considered. 

As can be seen from Table 1, a definitive conclusion 
based on the correlation coefficient would be difficult. 
Considering the heach totals, all four wave theories have 
comparable correlations; however, the regression equation 
for hyperbolic theory has the best slope and intercept com- 
bination, as further shown in Figure 5. 

There appears to be a wide range of scatter associated 
with the Santa Barbara data for which there is no simple 
explanation. Ecth, Torrey Pines and Leadbetter Beach, data 
sets were similarily collected and analyzed. While it is 
true that wave refractrcn diagrams were required for the 
Santa Barbera data, there is no reason to believe that there 
were significant errors in the computations. The offshore 
sensors from leadbetter Beach displayed littie scatter (Fig- 
ures 2-5) in the data; since these sensors were subject tc 
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the largest angles of incidence as veil as the largest cel- 
erities. There must be some other mechanism to explain the 
widespread scatter at the intermediate wave speeds. 

Guza and Thornton (19S2) have suggested that the pres- 
ence cf surf teat may be responsible for the differences 
between measured and theoretical celerities. Surf beat is a 
long reflected wave with a period cf the order cf several 
minutes, wavelength of the order of the surf zone width, and 
a maximum amplitude cf the order cf the swash at the beach 
face, decreasing offshore. Surf beat causes variations ir. 
the water depths, particular at the antinodes of the 
reflected wave, as perceived by the shorter sea-swell waves. 
This might explain seme cf the scatter or the data at Lead- 
better since the phase speed is depth dependent and sensor 
placement could have coincided with antinodal activity which 
would result in a lew frequency modulation of the water 
iepth above the current meters. Surf beat also induces ct- 



offshore 


velccit 


ies. 


particular at the 


nedes of 


re fleeted 


wa ves , 


whi ch 


can also affect the 


phase spe 



the sea-swell waves. Sensors placed near the nodes cf ler.g 
waves could experience large amplitude excursions in the 
waves’ horizontal velocities. These periodic differences in 
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the leng wave induced velocities could result in a dcp 
shifting cf the phase spectra and could aecount for sea 
in Lsadbetter data. 



IV 



. CONCL USIONS 

hyperbolic theory, which is an assymptotic fora of cnci- 
dal theory in shallow water, appears to give best agreement 
when all three data sets are collectively analyzed ana com- 
pared. LeMehaute et al (1968) similarily concluded cr.oidal 
theory gave the best comparison with laboratory data , while 
Dean (1965) found croidai theory is valid particularly in 
deep water hut is lacking for shallow water waves. Linear 
theory, on the ether hand, yields almost as good a compari- 
son and is computationally easier. The data did not allow 
for s proper test of wave theories in the deeper water and 
larger waveheight regime. 

Here and solitary wave theories generally cverpr edicted 
the measured celerities while hyperbolic and linear theory 
tended to slightly underpredict. The differences are 
attributable to the relative amplitude being dispersive, as 
predicted by various theories. Eore and solitary theory are 
strongly amplitude dispersive, hyperbolic theory is weakly 
dispersive, and linear theory has no amplitude dispersion. 
Therefore, it is concluded that the waves are best categor- 
ized as weakly amplitude dispersive. 
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Figure 1. Recicr. s of Aoplicability 

Torre y Pinas (+) Seven Mile Beach (o) 

L e edbe” ar ( a ) 
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MEASURED PHASE SPEED (CM/SEC) 

50 . 0150.0 250.0 350.0 450.0 550.0 650.0 750.0 350.0 950.0 1050.0 




Figure 2. Measured Wave Sceed vs. Linear Theory 

Tcrrey Pines (+J Seven Mile 3eacn (c) 

Leaabetner { a ) 
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MEASURED PHASE SPEED (CM/SEC) 

50.0150.0 250.0 350.0 450.0 550.0 650.0 750.0 850.0 950.0 1050.0 




SOLITRRY THEORY (CM/SEC) 



Figure 3. Measured Wave Sneed vs. Solidary Theory 
Tcrrsy Pines ( + f Seven Mile Beach (c) 

L eadberte r ( a ) 
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MEASURED PHASE SPEED (CM/SEC) 

50 . 0150.0 250.0 350.0 450.0 550.0 650.0 750.0 850.0 950.0 1050.0 




BORE THEORY (CM/SEC) 



Figure 4. Measured Wave Speed vs. 3ors Theory 

Torrsy Pines (+) Seven Mile 3each (o) 

L eadbe tte r { a ) 
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MEASURED PHASE SPEED (CM/SEC) 

50.0150.0 250.0 350.0 450.0 550.0 650.0 750.0 1350.0 950.0 1050.0 




HYPERBOLIC THEORY (CM/SEC) 

Fiaurs 5. Measured Wave Speed vs. Hyperbolic Thscr 
Torrey Pines (+f Seven Mile 3each (o 

L eadbe tre r ( a ) 
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TABLE I 

Statisorcal Comparison of Wave Theories 





LINEAR 


SOL IT A RY 


SORE 


HYPEREOLIC 


To rre y 
Pines 


Y=6< .3+.884X 


Y=28. 6+ . 855X 


Y=0 . 10 4 + . 847 X 


Y = 39 .0 + .9 19X 


Beach 


Ccr r=. 932 


Ccrr=.95 1 


Ccr r=. 955 


Cor r = . 9 20 


Se van 
Mile 


Y=62. 5+.759X 


Y = 36. 3+. 7 1 7x 


Y=-3 8. 0 + . 844 X 


Y= 57 . 0 + . b 67X 


Beach 


Corr=. 749 


Ccrr= .768 


Cor r = . 625 


Corr = . 6 97 


Le ad- 
be iter 


Y=74 . 5 + .89 6 X 


Y=30. 6+ . 862X 


Y=- 18. 9+ . 890 X 


Y = 5 1 .2+ .930X 


Be acn 


Ccr i: = . 897 


Ccrr= .89 1 


Cor r=. 883 


Ccrr=. 900 


Total 
of all 


Y=5 1.9 + .932X 


Y = 1 4. 5+. 886X 


Y=- 39.6+. 922X 


Y= 29 . 2+ . 9 75X 


Bsachs 


Ccr:: = . 9 37 


Ccrr= . 93 6 


Cor r= . 93 1 


Ccrr = . 9 3o 
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